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ABSTRACT 

We study the host galaxy properties of the tidal disruption object. Swift J164449.3-1-573451 using 
long-term optical to near-infrared (NIR) data. First, we decompose the galaxy surface brightness 
distribution and analyze the morphology of the host galaxy using high resolution HST WFC3 images. 

We conclude that the host galaxy is a bulge-dominant galaxy that is well described by a single 
Sersic model with Sersic index n = 3.43 ± 0.05. Adding a disk component, the bulge to total host 
galaxy flux ratio (B/T) is 0.83 ± 0.03, which still indicates a bulge-dominant galaxy. Second, we 
estimate multi-band fluxes of the host galaxy through long-term light curves. Our long-term NIR 
light curves reveal the pure host galaxy fluxes ~ 500 days after the burst. We fit spectral energy 
distribution (SED) models to the multi-band fluxes from the optical to NIR of the host galaxy and 
determine its properties. The stellar mass, the star formation rate, and the age of stellar population 
are log(M*/M 0 ) = 9.14lg;®Q, 0.03to'o3*^-63-o’43 Finally, we estimate the mass of 
the central super massive black hole which is responsible for the tidal disruption event. The black 
hole mass is estimated to be from Mbh “ AI*_bui and Mbh “ Tbui relations for the 

K band, although a smaller value of ^ 10^ Mq cannot be excluded convincingly if the host galaxy 
harbors a pseudobulge. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: photometry — galaxies: structure 
— techniques: photometric 


1. INTRODUCTION 

Swift J164449.3-f573451 (hereafter, Swift J1644-h57) 
was first discovered by the Swift Burst Alert Telescope 
(BATl at 12:57:45 UT o n 28 March 2011 (|Burrows et all 
I 2011 I iLevan et al.ll 2011 ^ . Some evidences suggest that 
Swift J1644-I-57 is a tidal disruption of a star by a su- 
permassive black hole (SMBH). This phenomenon trig¬ 
gered the BAT three t imes after the initial trigger dur¬ 
ing t he first few days ([Burrows et al.l[20Tlt ILevan et al.l 
1201 Ih . The late-time X-ray light curve was extended to 
a longer period by following the expected power-law de- 
cay for th e tidal disruption of a star, i.e., (e.g., 

lReeslll98^ . Finally, the source of X-ray, IR, and radio 
emissions were well matched up wit h the center of the 
host galaxy where a SMBH resides ([Levan et al.l 120111 : 

I Zanderer et al.ll2011h . 
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There have been many studies performed to under¬ 
stand the nature of this event, such as the characteristics 
of the star that was disrupted. Such a question is closely 
connected to the SMBH mass (Mbh)- The disrupti on of 
a solar-type star is possible for all Mbh < 10 ^ Mq di eei 
ligSSlICannizzo et al.lll990l : IBIoom et al.ll20Ilt . but com¬ 
pact stars like a white dwarf can be disrupted only if 
Mbh < 10® Mq (iKrolik &: Piranl[2011f) . If so, then this 
kind of event provides the interesting possibility of dis¬ 
covering intermediate-mass black holes. 

Unfortunately, there h as been contro v ersy c oncerning 
the mass of the SMBH. iBurrows et al.l ([20111 ) provided 


a rough estimate of the SMBH mass of 2 x 10^ A/q 
using a black hole mass - luminosity relation and the 
lower limit of ^ 10® Mp) b ased on the X-ray variabil¬ 
ity. Similarly, ILevan et al.l (1201111 estimated Mbh to be 
2 X 10® - 10^ Mq, derived from Al-band luminosity, but 
at that time, Al-band luminos ity contained a significant 
amount of the transient light. iMiller fc Giiltekinl (120111 ) 
utilized a relation between the black hole mass, the ra¬ 
dio luminosity, and the X-ray luminosity, and found Mbh 

10® ® M(Tt . Using a quasi-peri o dic os cillation resonance 
hypothesis. lAbramowicz fc Liul ( 2012f) provided an Mbh 
estimate of ~ 10® Mq. IKrolik fc PiranI ([201 Ih concluded 
that a white dwarf was tidally disrupted and the mass of 
SMBH is less than 10® Mq in light of the short timescales 
of the X-ray light curve. In summary, the Mbh esti¬ 
mates have centered around the two discrepant values 
of IO^’Mq and 10® Mq or less. A better understanding 
of the host galaxy properties is needed to clear up the 
situation. 

In order to more accurately estimate the SMBH mass 
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TABLE 1 

HST WFC3 Data Log 


Observation date (UT) 

MJD'* 

Days since trigger 

band 

Exptime (s) 

Magnitude (AB) 

2011-04-04 

55655.147614 

6.6 

F606W 

1260 

22.69±0.01 

2011-08-04 

55777.276876 

129 

F606W 

4160 

22.76±0.01 

2011-12-02 

55897.684390 

249 

F606W 

1113 

22.77±0.01 

2013-04-12 

56394.429204 

746 

F606W 

2600 

22.74±0.01 

2011-04-04 

55655.132654 

6.6 

F160W 

997 

20.68±0.01 

2011-08-04 

55777.257148 

129 

F160W 

1412 

21.09±0.01 

2011-12-02 

55897.702220 

249 

F160W 

1209 

21.22±0.02 

2013-04-12 

56394.295795 

746 

F160W 

2812 

21.55±0.02 

^ Exposure start time in Modified Julian Date (MJD) 






TABLE 2 





Spitzer IRAC Data Log 



Observation date (UT) 

MJD^ 

Days since trigger 

band 

Exptime (s) 

Magnitude (AB) 

2011-04-28 

55679.975316 

31.4 

3.6/im 

1250 

19.50±0.02 

2011-10-31 

55865.023052 

216 

3.6//m 

1253 

21.49±0.06 

2012-02-24 

55981.541644 

333 

3.6//m 

1252 

21.70±0.07 

2011-04-28 

55679.975316 

31.4 

4.5/im 

1250 

19.30±0.01 

2011-10-31 

55865.023052 

216 

4.5/im 

1253 

21.25±0.05 

2012-02-24 

55981.541644 

333 

4.5^m 

1252 

21.57±0.08 


^ MJD in UTC at data collection event (DCE) start 


and better constrain the properties of the host galaxy, we 
analyze the morphology and the surface brightness pro¬ 
file of the host galxaxy based on high-resolution Hubble 
Space Telescope {HST) images and estimate the multi¬ 
band fluxes of the host galaxy using our long-term mon¬ 
itoring data lasting more than 2.4 years. We fit the 
multi-band spectral energy distribution (SED) of the 
host galaxy luminosity with stellar population synthe¬ 
sis models, and then obtain the properties of the galaxy. 
Finally, we provide our best estimate of Mbh based on 
the host galaxy properties. 

This is the second of a series of two papers. In the 
first paper (M. Im et al. 2015 in preparation, hereafter, 
Iml5), we present the dataset of the long-term monitor¬ 
ing campaign and an analysis of the late-time light curve. 

Throughout the paper, we selected Hq — 70km 
s“^Mpc“^, r^A = 0.7, and V,m = 0.3 as cosmological 
parameters and adopt the AB magnitude system. 


2. OBSERVATIONS AND DATA 

We observed Swift J1644-I-57 using Wide Field Cam¬ 
era (WFCAM) on United Kingdom Infrared Telescope 
(UKIRT) for nearly 2.4 years following the burst as a 
part of our gamm a-ray burst (GR B) and transient ob¬ 
servation program (|Lee et al.ll2010ll . We observed inten¬ 
sively in the K band among Z, Y, J, H, and K bands of 
WFCAM. The number of epochs of K band data used 
for thie analysis is 101 and the last data were observed 
at 884.7 days after the initial BAT trigger. The num¬ 
bers of epochs for the Y, J, 7J-band data are 3,15, 28 and 
the last data were observed at At = 712.1, 710.1, and 


884.7 days, respectively, where At is the number of days 
since the initial BAT trigger. We have only one epoch 
of data for the UKIRT Z band which was observed at 
At = 723.0 days. 

We also observed Swift J1644-I-57 usin g Camera 
for QUasars in EArly uNiverse (CQUEA N; iKim et HI 
120111 iPark et ^ 120121 iLim et all 1201311 on the 2.1m 
Otto-Struve telescope of the McDonald Observatory in 
g,r,i,z, and Y bands. The numbers of epochs for the 
g, r, i, z, and Y band data are 2, 2,14,14, and 2 and the 
last data were observed at At = 25.7,217.5,526.6,526.6, 
and 25.8 days respectively. The UKIRT and CQUEAN 
observation logs and photometry results are described in 
Iml5. 

In a ddit ion, we also used data from iBurrows et al.l 
(|2011ll and I Levan et al.l (1201111 for the earlier optical and 
near-infrared (NIR) data. 

Morphology analysis requires high-resolution images 
because this object is so compact that it is virtually 
a point source in the UKIRT and CQUEAN images. 
For the high-resolution imageSjWe obtained HST WFC3 
multi-drizzled, stacked imaged available in the MAST 
database. We used F606W-, F160W-band data and the 
number of epochs in each two band is four. These data 
were observed at At = 6.6,129, 249, and 746 days. The 
HST WFC3 data are summarized in Table [TJ 

To supplement the NIR observation data, we used 

® Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space Tele¬ 
scope Science Institute. STScI is operated by the Association of 
Universities for Research in Astronomy, Inc. under NASA contract 
NAS 5-26555. 
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Fig. 1. — Images of the host galaxy in F606W, the best-fit two-dimensional models from GALFIT, and the residuals for the single Sersic 
component model and the Sersic bulge -|- exponential disk model. Right panels show one-dimensional profiles of the host galaxy and each 
model component. 


the Spitzer IRAC 3.6, 4.5/.tm post basic calibrated data 
(PBCD) from the NASA/IPAC Infrared Science Archive. 
These were observed at At = 31.4, 216.5, and 333.0 days. 
A log of the Spitzer IRAC 3.6, 4.5/rm data are shown in 
Tabled 

The flux me asurements were perform ed by SExtrac- 
tor software^ (iBertin fc Arnouti 1199611 except f or the 
HST images for which GALFIT (jPeng et al.ll2010f) mod¬ 
els were used for the flux measurements. 

For the X-ray data, we used Swift/XRT data taken 

® We used aperture magnitudes with aperture correction. 


from the Swift archive and XMM-Newton datc(3 from 
the XMM-Newton Science Archive. 

3. MORPHOLOGY OF THE HOST GALAXY 

We analyzed the surface brightness profile of the host 
galaxy, in order to determine the bulge fraction and its 
nature. We used the HST images and GALFIT software 
to fit two-dimensional models to the light distribution of 

Based on observations obtained with XMM-Newton, an ESA 
science mission with instruments and contributions directly funded 
by ESA Member States and NASA 












































4 


Yoon et al. 


Oriainal 5 kpc 

1" 

• 

Model 

• 

Residual 

4 

Exponential Disk 



0.4 0.6 0.8 1.0 1.2 

(Semi major axis)^''' [arcsec^^''] 


Oriainal 5 kpc 

1" 

Model 

A 

A 

• 

9 

Residual 



Sersic Bulge (n=1) 

+ Exponential Disk 



0.4 0.6 0.8 1.0 1.2 

(Semi major axis)^^'' [arcsec^^''] 


Fig. 2.— Same as Figure [T] but for the single exponential disk model and the double exponential profile model. 


the host galaxy. To construct the point spread function 
(PSF), we selected ^ 5 isolated stars with signal-to-noise 
ratios > 300 in the vicinity of Swift J1644+57, and co¬ 
added them. 

We used error images that are created by GALFIT 
for the fitting. For GALFIT to create the error image 
properly, we modified the unit of ADU and the image 
header values such that GAIN x ADU x NCOMBINE = 
[electrons] as recommended in GALFIT websitJ*^. 

A crucial factor affecting the fitting results is back¬ 
ground subtraction. For the background determination, 
we set 6 annuli with the radii logarithmically increasing 


between 2.5 and 9 times the radius of an ellipse for which 
pixel values are I.5 ct of the background noise. We cen¬ 
tered the annuli on the center of the host galaxy, set the 
minimum width of the annuli to be ~ 1.3 arcsec (^ 33 
pixel) for F606W images and ^ 2 arcsec (^ 16 pixel) for 
F160W images, and augmented the widths in step with 
the logarithmically growing radii. We then derived the 
mean pixel values of each annulus. Finally, we adopted 
their mean value as the background value. 

Our surface brightness fit was carried out using a deep, 
stacked image of the data taken with F606W at At = 
129, 249, and 749 days. It has been known that the 


http: //users.obs. camegiescience.edu/peng/work/galfit/TOP10.html 


transient component is negligible in the optical bands 
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Fig. 3. — Flux fractions of the GALFIT models as a function of 
time. We used the model which consists of a single Sersic bulge 
with n = 3.43 and a point source component. The upper panel 
shows the results for the F606W-band images, while the lower panel 
shows the results for the F160W-band images. Magnitudes of the 
model components are also shown. 

bluer than i even at the early time (|Burrows et alJ[2011l : 
iLevan et alJl20ljT) . The use of the stacked, late-time im¬ 
age in the F606W band makes the transient component 
more negligible. On the other hand, NIR-bands, includ¬ 
ing F160W (similar to H band of WFCAM), are known 
to contain a significant transient component which may 
affect the host galaxy analysis. Furthermore, the spatial 
resolution of the F606W images is better by a factor of 
three than that of F160W, which greatly helps the sur¬ 
face brightness fitting. The other HST data were also 
analyzed to understand the importance of the transient 
component, and the results for the transient component 

are presented later in this section. __ 

For the galaxy models, we used the S ersic (ISersidllQ^ . 
de Vaucouleurs (|de Vaucouleursl[l948n . and exponential 
disk profiles or a combination of thereof. The Sersic pro¬ 
file is described as 


S(r) 


Se^Xp 


— K 



where Eg is the surface brightness at the effective radius 
Tg, and n is the Sersic index, k is a variable parameter 


denpendent on n, where n = 4 and 1 correspond to the de 
Vaucouleurs and exponential profiles, respectively. Al¬ 
though n = 4 is commonly quoted for the ellipticals and 
classical bulges, the Sersic index of ellipticals and classi¬ 
cal bulges can assume a value in the range 2 < n < 6 , 
wher eas pseudobulges have n <2 (jFisher fc Drorvl[2nn8l 

IMnli . 

All of the model parameters such as ellipticity and cen¬ 
ter positions of the different components, were set free in 
the fitting procedure. 

Figures [T] and [2] show images of the host galaxy, 
the two-dimensional models, and the residuals (i.e., the 
model subtracted images), for four different models: (1) 
a single Sersic; (2) a S&sic bulge -I- exponential disk; (3) 
an exponential disk; and (4) a double exponential profile 
models. The figures also show one-dimensional profiles 
(along the major axis) of the host galaxy and those of 
each model component, which are converted through the 
IRAEI3 ellipse task. In addition to the profiles, the 
differences between the data and the model profiles are 
shown. The results of each fit are summarized in Ta¬ 
ble [3l Both the single Sersic model with n = 3.43 ± 0.05 
and the Sersic bulge with n = 3.39 ±0.11 ± exponential 
disk model fit the data well {xv 1-2 - 1.3). When the 
disk component is added, the bulge to total host galaxy 
flux ratio (B/T) is 0.83 ± 0.03. On the other hand, the 
single exponential disk model provides a poor fit to the 
data as shown in Figure 2 and with Xu = 6.54. The 
double exponential profile model fits the data nearly as 
well as the single Sersic model and the Sersic bulge±disk 
model in terms of xt- However, the analysis of the 
one-dimensional surface brightness profile shows that the 
model does not follow the outer part of the profile well, 
demonstrating a relatively steeper decline than that of 
the single Sersic model and the Sersic bulge ± exponen¬ 
tial disk model. This model gives B/T = 0.36, suggest¬ 
ing a significant bulge component. Therefore, we con¬ 
clude that the host galaxy of Swift J1644±57 is bulge- 
dominant. We also conclude that the bulge is likely to 
have a Sersic index higher than 3 regardless of the exis¬ 
tence of the disk. Thi s value corresponds to the range 
of the classical bulges (iFisher fc Drorvl[2008L 120101) . We 
cannot completely exclude the case where the bulge is 
pseudobulge with n ~ 1, but even in this case, the ob¬ 
ject has a significant bulge. 

Additionally we fit the observed surface brightness pro¬ 
file with a single de Vaucouleurs bulge model and a de 
Vaucouleurs bulge ± exponential disk model. The results 
of these fits are nearly identical to that of the single Sersic 
and the Sersic bulge±disk models. 

To estimate the transient component flux, we fit all 
the F606W and F160W images with a model containing 
both the point source (transient) and the host galaxy 
components. Here, we adopt a single Sersic profile with 
a fixed Sersic index (n = 3.43) for the host galaxy com¬ 
ponent, and a PSF profile for the transient component. 
The compactness of the host galaxy and the bright tran¬ 
sient component in the F160W images create a serious 
degeneracy, particularly between the effective radius 

IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under a cooperative agreement 
with the National Science Foundation. 
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TABLE 3 

Surface Brightness Fitting Result 


Galaxy model 
(1) 


Bulge 


Disk 

B/T 

( 7 ) 

mt [AB] 

(8) 

xl 

(9) 

mb [AB] 

(2) 

n 

( 3 ) 

Teff [kpc] 

( 4 ) 

md [AB] 

( 5 ) 

rs[kpc] 

(6) 

B 

22.77 ±0.01 

3.43 ± 0.05 

1.01 ±0.01 




22.77 ±0.01 

1.322 

B+D 

23.01 ±0.03 

3.39 ±0.11 

0.79 ± 0.03 

24.75 ±0.14 

1.21 ± 0.07 

0.83 ±0.03 

22.81 ±0.03 

1.223 

B(n = 4) 

22.72 ±0.00 

4 (fixed) 

1.09 ±0.01 




22.72 ±0.00 

1.275 

D 




23.07 ±0.00 

0.47 ±0.00 


23.07 ±0.00 

6.547 

B(n = 4)+D 

23.01 ±0.03 

4 (fixed) 

0.85 ±0.02 

24.64 ±0.10 

0.95 ± 0.03 

0.82 ±0.03 

22.79 ±0.03 

1.228 

S(n = 1)+D 

23.99 ±0.01 

1 (fixed) 

0.30 ±0.00 

23.35 ±0.00 

0.92 ± 0.01 

0.36 ±0.00 

22.87 ±0.00 

1.351 


Note. — Column 1: galaxy model for the two-dimensional fitting. B: Sersic bulge, D: exponential disk , B(n = 4): de Vaucouleurs 
bulge. S(n = 1): Sersic profile with fixed n — 1 (exponential profile). Column 2: AB magnitude of the bulge component. Column 3: 
Sersic index for the bulge model. Column 4: effective radius. Column 5: AB magnitude of the disk component. Column 6: scale length 
of the disk component. Column 7: bulge to total light ratio. Column 8: total magnitude. Column 9: reduced foi" th® fitting model 
defined as 

2 _ 1 {fdaLta-i^j v) /model v)) 

CT(a:,y)2 

where fdatai^-, v) a^nd /model(^» v) mean the input data and the model images, respectively. A^dof is the degree of freedom. cr(a:, y) is the 
error image. Here, sum is only over all nx and ny pixels satisfying 1.5cr of the background noise. 
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Fig. 4.— N I R and the X-ray l i ght cu rves of the Swift J1644-I-57. The light curves of Y,J,H,K bands with the early data from 
IBurrows et al.l II2Q1U) . [Levan et al.l (120111 ). and Swift/X.KT 0.3 - lOkeV are shown. All of the light curves have similar shapes except that 
the X-ray light curve is ~15 days ahead of the NIR light curves. After ^^500 days the X-ray emission was rapidly declined as shown in with 
star mark for the last X-ray data from the Swift/X.Kr. 

and the Sersic index, when fitting multi-component band. The flux fractions of the models as a function of 

models. To alleviate the degeneracy, we fixed the Sersic time are shown in Figure [31 In the case of F606W, the 

index to be n = 3.43, similar to that of the F606W flux fraction from the transient component is very small 
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Fig. 5.— Swift and XMM-Newton data as well as recent 
Chandra observation of Swift J1644+57. The X-ray flux abruptly 
declined after ~500 days since the BAT trigger. 



Days since BAT trigger 

Fig. 6.— H- and iT-band light curves in linear scale. Some of 
the data points at very early times are cut in order to highlight the 
late-time light curves. The magnitudes of the latest H, K bands 
converge to single values, suggesting that the transient component 
has disappeared at At > 500 days. 

or nonexistent, while the transient component is very 
bright in the earliest F160W band, even brighter than 
the entire host galaxy. The result reflects a very red 
color to the transient well, and justifies the exclusion 
of the point source component in the late-time F606W 
images during the host galaxy analysis. The point 
source contribution declines rapidly as time goes on in 
F160W, but it contributes to the total flux of the object 
until around At = 750 days. On the other hand, the 
fluxes of the host galaxy component are almost constant 
in both bands over the entire period. The magnitude of 
host galaxy in the F160W band is ~ 21.75 mag, and as 
we shall see in the next section, this is the same as for 


the last data point of the UKIRT 77-band light curve, 
suggesting that the flux of the last data point in the 
NIR light curve represents the host galaxy flux. 

4. LIGHT CURVES 

In this section, we show long-term observation results 
and estimate multi-band fluxes of the host galaxy of 
Swift J1644-I-57. Figured shows the Y, J, H, and K light 
curves. The gray data points in the background show 
Swift/XRT 0.3 - lOkeV data. The J,H,K light curves 
resemble each other. The NIR light curves rapidly de¬ 
cline until At ~ 10 days, turn up again with a second 
peak at At ~ 30 days, and decline again steadily. The 
behaviors of these NIR light curves are very similar to 
that of the X-ray light curve except that the X-ray light 
curve appears shifted ahead of the NIR light curves at 
a time of ^15 days. The similar shapes of these light 
curves indicate that the origins of the X-ray emission 
and NIR emission are related to each other. On the other 
hand, the time gap between these two emissons denotes 
that the X-ray source and NIR sourc e are separated from 
each other as much as the time gap. iBloom et al.l (j2011f ) 
suggested that X-ray source is in the close vicinity of the 
black hole due to the fact that the X-ray emission shows 
very rapid, high variability, while the IR and the radio 
sources are located a large distance from the black hole 
on account of the relatively smooth and small variability. 
They argued that the jet generated by black hole collides 
with the surrounding medium where the electrons are ac¬ 
celerated by the jet. These high-speed electrons emit the 
IR to radio photons through synchrotron radiation. 

The jet seems to be nearly turned off at At ~ 500 
days in light of the fact that there is an abrupt decrease 
in flux of a factor of ^ 10 or more, which can be seen 
in all the Swift/XRT, Chandra, and XMM-N ewton data 
(| Levan fc Taiivii]r201^ iZauderer et all 1201 3D . This late- 
stage turn-off of X-ray flux is also shown in Figure [SJ 
If the jet was turned off, then the transient components 
of the NIR fluxes must be quenched following the X-ray 
flux, and it is expected that the fluxes of the pure host 
galaxy of Swift J1644-I-57 were revealed at that time. As 
we can see in Figure IHl which shows the 77-, K- bands 
light curves in linear scale, the fluxes of the 77, K bands 
converge to single values at late-time. Furthermore, the 
latest 77-band magnitude is nearly the same as that of 
the host galaxy of the F160W-band images, shown in the 
results of the model fitting in ij3l This evidence indicates 
that it is reasonable to regard the NIR {Y, J, 77, and K 
band) fluxes of the last data, taken at At =~ 700 or 884 
days, as those of the pure host galaxy. 

The left panel of Figure [7] shows the light curves of 
the CQUEAN i- and z-band and the UKIRT Z-band 
data. In the case of the z band, the fluxes from the 
object slightly decrease with time. We also regard the 
flux of the last data, that is the UKIRT Z-band data, as 
the Z-band flux of the host galaxy since it is observed 
far beyond expected quenching time of jet. On the other 
hand, i-band fluxes are virtually constant, demonstrating 
that the transient components are basically non-existent 
in the i or bluer bands (Figure 151 iBurrows et al.l [20111 : 
iLevan et al.ll20llh . We take the flux of the last data 
of the i band as that of the host galaxy. The number 
of CQUEAN g- and r-band data points are scarce com- 
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Fig. 7.— Light curves of CQUEAN i and 2 : bands, UKIRT Z-band data, and Spitzer IRAC 3.6/i.m and 4.5/im bands. Note that the 
2 :-band flux decreases with time, while the i-band flux is almost constant with time. Considerable magnitude changes in the Spitzer IRAC 
3.6/i.m and 4.5/i.m bands can be seen in the right plot. 



Fig. 8.— Temporal change of the SED of Swift J1644+57. The 
fluxes in the redder bands show the substantial changes with time, 
whereas the bluer-band fluxes do not vary with time. 

pared to the NIR data and the last data were observed 
at early time {At = 25.7, 217.5 days, respectively). How¬ 
ever, there is little or no change between the very early- 
time magnitudes from ILevan et al.l (1201111 and our g- and 
r-band magnitudes in the same way as the fluxes of the 
i and F606W bands. Therefore, we consider the g- and 
r-band fluxes of the last epoch data as those from the 
host galaxy. Furthermore, we also consider the magni¬ 
tude of the single Sersic model of the stacked HST WFC3 
F606W image as that of the host galaxy since the point 
source contribution to whole flux is negligible. 

The right panel of Figure 0 shows the light curves of 
the Spitzer IRAC 3.6/rm and 4.5/rm bands. The magni¬ 
tude changes in these bands are the more significant than 


TABLE 4 

Magnitudes of Host Galaxy 


Band 

Magnitude (AB) 

B 

24.24-1-0.10 ILevan et al. 20111 

9 

23.67±0.19 

F606W 

22.72±0.01 

r 

22.73±0.06 

i 

22.25±0.04 

z 

22.16±0.06 

Y 

22.18±0.07 

J 

21.96±0.08 

H 

21.74±0.12 

K 

21.55±0.10 

3.6/.im 

21.70±0.07 (Including transient) 

4.5pm 

21.57it0.08 (Including transient) 

Note. 

— Magnitudes are cor- 

rected by 

Galactic extinction based on 

ISchlaflv Finkbeiner] l|201lh. 


those for the other optical/NIR bands. The fluxes from 
the transient component seem to be non-negligible even 
in the last epoch data {At = 333 days), because the last 
IRAC epochs were still in the rapidly decreasing phase. 
Therefore, we consider the fluxes of the last epoch IRAC 
data to be the upper limit fluxes of the host galaxy. 

The multi-band magnitudes of the host galaxy of Swift 
J1644-I-57 are shown in Table IH We took the Galac¬ 
tic extm£don_jrfj3hotometricdata into account based 
on iSchlafiv fc Finkbeiireij (|2011ll. We add ed the H-band 
photometric data from lLevan et aP (|2011h to expand the 
data points to the short wavelength band. 

The temporal change of the observed SED of Swift 
J1644-I-57 are summarized in Figure [5] The fluxes 
in redder bands show substantial changes with time. 
Meanwhile, bluer band fluxes are constant. This red 
feature of the transient has been suggested to be due to 
dust extinction. From previous studies, it is known that 
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TABLE 5 

Input Parameters for SED Fitting 


Parameter 

Value 

r (e-folding time scale of stellar population) 
t (age of stellar population) 
IMF (initial mass function) 

Z (metallicity) 
Extinction law 

Ay (V-band attenuation for stellar population in magnitude) 

6.5 < log[T/yr] < 11.0 with a step size of 0.1 
8.0 < log[t/yr] < 10.3 with a step size of 0.1 
Saloeter ('1955) 

0.004, 0.008, 0.020, 0.050 

Calzetti et al. ('2000) 

0.0 < Ay < 3.0 with a step size of 0.1 


TABLE 6 

Best-fit Parameters from SED Fitting 


Parameter Value 

Stellar mass [log(M*/M 0 )] 9.14^q'jq 

SFR [Mq/yt] 0.03l“;g« 

Specific SFR [log(sSFR/yr-l)] -10.62t^®° 
t [Gyr] 0.63l“;|5 

r [Gyr] 

Au 0.00t“:^J 


Z 0.050t°;°°° 
1,64 


X. / iim [rest-frame] 

0.5 1.0 3.0 



Fig. 9.— SED fitting result of the host galaxy of Swift 91644+57. 
Plotted are the multi-band fluxes of the host galaxy in Table [4| 
Two Spitzer data are the upper limit fluxes. The value for the 
fit is 1.64. 

the hydrogen column density of the source of X-ray is 
large (Vh ^ lO^^cm”^), meaning the line of sight to the 


SMBH has a very larg 

e extinction value (Ay = 4.5 - 

10) ( 

Bloom et al. 

2011| 

Burrows et al.ll201 It ILevan et al.l 


: IShao et al.l 

2011 

Saxton et al.ll2012fl. 


5. SED EITTING 

We performed SED model htting of the multi-band 
fluxes of host galaxy of Swift J1644-1-57 in order to de¬ 
termine the properties of the host galaxy such as the stel¬ 
lar mass (M*) which is an important parameter for the 
A4 bh estimation. We utilized the code Fitting and As¬ 


sessment of Synthetic Templates (FAST F^ iKriek et al.l 
l2009f) . which is a public SED fitting tool for the inves¬ 
tigation of the galaxy properties using the photometric 
data ranging from UV to IR. The code is based on the 
IDL and fits of UV to IR stellar population templates to 
photometric data or galaxy spectra. FAST runs with the 
method of fitting and using stellar population grids 
to derive the best-fit model and its parameters. 

We used the 2003 version of the Bru zual & Chariot 
(BC03) model (jBruzual fc Charlotl l2003fl for the stellar 
population model. There are thr ee initial mass f unctions 
(IMFs) available in the FAST (jSalpeteil 119551 : iKroupal 
1200 H l^abried [2003H . We chose the Salpeter IMF. To 
define the star formation history (SFH), we assumed an 
exponentially decreasing star formation rate (SFR). The 
stellar population was modeled with the e-folding time 
scales, 6.5 < log[r/yr] < 11.0 with a step size of 0.1 and 
ages of 8.0 < log[</yr] < 10.3 with a step size of 0.1. 
We used several metallidty values such as Z = 0.004, 
0.008, 0.02, and 0.05. The model SEDs were attenuated 
by dust, for which w e used attenuation curves based on 
iCalzetti et al.l (j200r)l] . We adopted 0.0 < Ay < 3.0 with 
a step size of 0.1. All the input parameters for the SED 
fitting are summarized in Table [5l 

Figure E] shows the SED fitting result. The two Spitzer 
data were treated as the upper limit fluxes of the host 
galaxy. The estimated stellar mass of the host galaxy is 
log(M*/M 0 ) = 9.141q ]q. The e-folding time scale is r = 
O.lOto.io Gyr and the age of the stellar population is t = 
0.63to,43 Gyr. The SFR of galaxy is 0 . 031 q;q 3 A/ 0 /yr, 
and the spe c ific SF R is log(sSFR/yr-i) = -10.621)^^°. 
iLevan et al.l (|2011ll derived SFR of 0.3 - 0.7 Mq/jt from 
the Ha and [O H] emission line luminosities. The value 
of 0.3 Mq/jt from Ha is consistent with our la upper 
limit. The SFR from [O H] line {0.7 Mq/jt) is about 
twice larger but the [O H] line based SFRs are known to 
be dependent on phys ical condition such a s the reddening 
and metallicity (e.g.. iKewlev et ani2004D . and less reli¬ 
able than Ha based SFRs. Another possible cause of the 
discrepancy is the different timescales that are probed by 
different SFR indicators ( emission line indica tors prob¬ 
ing recent star formation). ILevan et al.l (j 2011 f ) estimated 
E{B — V)gas = —0.01±0.15mag, i.e., no extinction using 
the intensity ratio of Ha and H/3 lines. This is consis¬ 
tent with our SED fitting result Ay = O.OOt^oo- The 
value for the fit is 1.64. 

The host galaxy of Swift J1644-I-57 is a low mass, low 
SFR galaxy with a low extinction. Also it seems to have 
experienced a rapid decline of SFR not v ery long ago. 
This fits in well with a recent suggestion bv lArcavi et al.l 

http://astro.berkeley.edu/~mariska/FAST.html 
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(|2014f) that host galaxies of tidal disruption events are 
E+A galaxies with < 1 Gyr stellar population and low 
or no SFRs. 

The output parameters are given in Table |6l The er¬ 
rors correspond to Ict confidence intervals derived from 
100 times Monte Carlo simulations in which the input 
photometric data are changed according to their errors. 

We also tried the Chabrier IMF instead of the Salpeter 
IMF for the fit. The change of the IMF influenced to the 
stellar mass, decreasing the stellar mass by ^ 0.25 dex. 

We also fitted the SED model with the IMarastonl 
(|2005f) stellar population instead of the BC03 model. We 
set the input parameter ranges identical to the case of 
the BC03 stellar populations. The results were nearly 
identical to the BC03 result. 

Our analysis of the host galaxy shows that the host 
galaxy is a bulge-dominated and nearly extinction 
free (Ay ~ 0 mag). On the other hand, the spectral 
properties of the nuclear transient suggests a high 
extinction (Ay 6 mag). These two facts may appear 
contradictory, but we note that a significant amount 
of dust can be found easily in nuclear region of bulge- 
dominated galaxies when their nuclei are acitve. For 
example, hosts of luminous AG Ns are mostly early -type, 
bulge-dominated galaxies (e.g.. iHong et alJ 1201511 . and 
such AGNs are known to contain a significant amount 
of dust in nuclear region as a form of hot or warm dusty 
torus (e.g., iKim et a].l[2(rill . 


Our^re^(MBH-LK,bui) 
Our result (Mbh-M| 3 ^,|) 
Our result (if pseudobulge) 


—■ Krolik & Piran (2011) 

-Abramowicz & Liu (2012) 


Millers Gultekin (2011) 

Levan et al. (2011) 


Burrows et al. (2011) 


4 5 6 7 8 

loglMen / M@] 


6 . DISCUSSION ON BLACK HOLE MASS 

Our results on the properties of the host galaxy of Swift 
J1644-I-57 can be summarized as follows. It is a bulge- 
dominant galaxy (B/T=0.83 ± 0.03). The mass of the 
host galaxy is somewhat low at 10®'^^ Mq, even though 
the galaxy is bulge-dominated. Now we estimate the 
mass of the SMBH that played the main role in the tran¬ 
sient phenomenon. 

It is now generally accepted that the SMBHs (10® - 
10^® M 0 ) reside in the bulges of all massive galaxies. 
Tight scaling relations have been derived between SMBH 
mass and several physical properties of the bulges (veloc¬ 
ity disper sion, mass, lumin o sity, e tc .) in many previous 


studies (IMaeorrian et al.l 119981: iFerrarese & Merritt 

200 c 

: IGebhardt et al.l 120001: McLure & Bunion 

2002 

; IMarconi & HuntI 120031: iHarinff & Rixl 

2004; 

Aller & Richstonel 120071: iHookins et all 

2007; 

Giiltekin et al. 

20091: Kormendv & Bended 

2009; 


Sani et al.l 120111 : iKormendv fc Hoi I2013D . Some argue 


that ellipticals and classical bulges follow identical rela¬ 
tions, while the pseudobulges follow a s o mewhat different 
relati o n with large sca tter (iHrj I 2 OO 8 I: iKormendv et al.l 
I 2 OIII : iSani et ahl 120111 : IKormendv fc Hoi 120131 ). We 
conclude from the best-fit galaxy models that the host 
galaxy of Swift J1644-1-57 has a classical bulge, and have 
also found a minor possibility of the pseudobulge with 
B/T=0.36. For now, we consider only the best model, 
that is, the case of the host galaxy having a classical 
bulge and being bulge-dominant. 

In order to esimate the central SMBH mass in the 
host galaxy of Swift J1644-I-57, we used the scaling re¬ 
lation between Mbh and the stellar mass of the bulge 
(Ar*,bui)- We expect that a large part of the s tellar mass 
derived in 1j5]belongs to the bulge component. ISani et ^ 


Fig. 10.— Results on the black hole mass from this work and 
previous studies. The red circles and dashed line represent our 
results, while the black squares and lines denote the results from 
previous studies. Error bars correspond to deviation of Icr and the 
arrow indicates the upper bound value. 


(|2011ll present the Mbh - M*_bui relation, where M*^bui 
is directly obtained from the bulge luminosity (Lbui) of 
Spitzer 3.6/rm and the calibrated M*^bui “ Tbui relation. 
They excluded pseudobulges when constructing the rela¬ 
tion. The relation is 


log(A/BH/M©) =a + f3x [log(M*_bui/MQ) - 11], (1) 

where a = 8.16 ±0.06, jS = 0.79 ±0.08, and the intrinsic 
scatter is 0.38±0.05 . The estimated mass of the SMBH is 
2 q6.7±o. 4 based on the stellar mass of the host galaxy 
and the above relation. If we consider the B/T=0.83 and 
assume that the mass-to-light ratio is constant in the 
bulge and disk, the stellar mass is decreased by ^ 0.1 
dex. It leads to a decrease in Mbh by ~ 0.1 dex. 

The tight scaling relations between Mbh and host 
galaxy properties suggest a close link between the 
SMBH growth and the galaxy evolution. There may 


which there have been various studies 

Treu et al. 

2004 

McLure et al.ll2006l: iShields et al.ll2006; 

Woo et al. 

2004 

Salviander et al. 20071: Treu et al.ll2007 

Woo et al. 

20081 

Jahnke et al. 20091: Bennert et al.l 20 

Ol Decarli et al.l 

20101:IMerloni et al.ll2010l: IBennert et al.ll2011l). The evo- 


lution of the scaling relations is still controversial in terms 
of the selection effects in the high redshift regime. Nev¬ 
ertheless, we can consider a case where the growth of 
Mbh happened ahead of the assembly of t he stellar mass 
as sug gested by some of these studies. iBennert et al.l 
(|201lh suggest the redshift evolution out to z ~ 2, in the 
form of Mbh/M^^buI oc (1 ± 2 ) 1 ®®=*=®-®®, based on 11 X- 
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ray-selected broadline AGNs. Gonsidering this evolution 
effect and the redshift of Swift J1644-I-57 z = 0.35, the 
mass of the SMBH could be larger by ~ 0.3 dex. 

We also estimated Mbh through th e Mrh - K band 
lumin osity of bulge (Lx^bui) relation in IKormendv fc li3 
(tMl . assuming that most of the NIR fluxes come 
from the bulge. Their relation is much improved com¬ 
pared with previous studies in view of serveral things. 
They excluded galaxies with black hole monsters which 
have over-massive SMBHs despite having relatively small 
bulges and ellipticals. They also omitted galaxies with 
black hole masses that are measured based on the kine¬ 
matics of ionized gas without taking line widths into 
account, since this method may yield underestimated 
masses. Galaxies in the process of merging generally have 
low-mass black holes for their luminosities. Thus they ex¬ 
cluded these galaxies from the relation. They also did not 
includ e pseudobulges. The relation of iKormendv fc H^ 
(IM^ is 

log(MBH/10'' Mq) =-a-13 X (Mk, bui + 24.21), (2) 

where the a is 0.265 ± 0.050, the [3 is 0.488 ± 0.033, and 
intrinsic scatters is 0.30. Mk, hui is the AT-band absolute 
magnitude of the bulge based on the photometric system 
of 2MASS. Using the best-fit SED derived earlier, we 
applied K-correction and evolution correction. We find 
Mk = —19.84 Vega mag including the evolutionary cor¬ 
rection of 1.38 mag, which is derived by the difference in 
the AT-band magnitude between 0.63 Gyr old population 
as in Table m and 4.5 Gyr old population, which is the 
age of this host galaxy in the local universe. Then, the 
mass of the SMBH is estimated to be 10 ® ®=*=°'^ Mq. If we 
take B/T=0.83 into account and assume that this value 
is also applicable to the NIR bands, then Mrh decreases 
by ^ 0.1 dex. If we use an Mk value that has not been 
corrected for evolution, then Mrh becomes ^ 0.7 dex 
larger. 

The tidal disruption of normal stars by a black hole 
is not possible for Mbh > 10® Mq, since the tidal ra¬ 
dius where the disrupt ion can occur is located inside the 
Schwarzschild rad ius (jReesI [19^ iCannizzo et al.l IT990I : 
iBloom et al.ll20lHl . Our results on the mass of the SMBH 
satisfy the condition for a tidal disruption event. 

Although we favor the model in which the host galaxy 
of Swift J1644-I-57 is a classical bulge, our analysis 
shows that it could be a galaxy with a pseudobulge with 
B/T=0.36 (Table [3|). If so, it is rather difficult to ob¬ 
tain an Mbh value, since the scaling relation is not well 
established for pseudobulges, especially in the low mass 
range of M* ~ 10® Mq for the host galaxy. Several works 
have shown that the Mbh ~ host galaxy scaling relations 
are weak or zero with a large scatter for pseudobulges. 
Over the M*_bui range of 10®'® to 10®® ® Mq where such 
a relation has been studied, Mbh cai i have any value 
between 10® to 10® Mq (Figure 21 of iKormendv fc: h 3 
1201311 . To reach down to M*^bui ~ 5 x 10® Mq as implied 
from the pseudobulge fit of our data, currently one can 
barely do so by re l ying on results from low mass AGNs 
(iBarth et al.ll2005l : iGreene et al.ll200§ : iJiang et 311120111 : 
iXiao et al.ll2011ll In such a case, an Mbh value between 
10® t o 10®-® Mq is possible (Figure 32 of iKormendv fc: H^ 
1201311 . Overall, if the host galaxy harbors a pseudob¬ 
ulge, then we can only loosely constrain Mbh to have 


a value between 10® to 10 ^’Mq considering our current 
poor knowledge of the Mbh value in pseudobulges. 

It is also known that a small fraction of pseudobulges 
have a Sersic index of n > 3. The best example is Pox 
52, for which n ^ 3.6 - 4.3, Mrh 2 x M^, an d 
M* ~ 10® Mq (jBarth et al.ll2004riThornton et al.ll2008ll . 
Therefore, even if we accept the Sersic index of n = 3.43 
as the best-fit result, we need to keep this kind of caveat 
in mind. 

Figure [TO] shows our overall results on Mrh and 
the results from the previous studies we mentioned in 
SI It shows that our favorite results are_comgatib]e 
with t he pr evious rough estim ates from iBurrows et al.l 
llml and I Levan et al.1 (1201111 . who also used scaling 
relations. However, our results are improved compared 
to the previous results, by revealing that the host 
galaxy has a significant bulge component through a 
two-dimensional bulge -I- disk decomposition of the 
surface brightness profile, and removing the transient 
component in NIR light using a long-term light curve. 
The Mbh limit could be much looser (the dashed line) 
if the host galaxy harbors a pseudobulge. A critical test 
of the pseudobulge model would be to obtain a deep, 
high-resolution image to see how the surface bright¬ 
ness profile behaves at the outer region of the host galaxy. 


7. SUMMARY 

We investigated the host galaxy properties of tidal dis¬ 
ruption event. Swift J1644-I-57 through morphology anal¬ 
ysis, light curve analysis, and SED fitting. We also es¬ 
timated Mbh which played the main role of this phe¬ 
nomenon, through scaling relations. 

We decomposed the surface brightness profile of the 
host galaxy based on high-resolution HST WFG3 images. 
We found that the host galaxy of Swift J1644-1-57 is a 
bulge-dominated galaxy which is well described by a sin¬ 
gle Sersic model with the Sersic index, n = 3.43±0.05. If 
we add a disk component, the bulge to total host galaxy 
flux ratio (B/T) is 0.83 ± 0.03, still indicating a bulge- 
dominant galaxy. We conclude that the host galaxy of 
Swift J1644-I-57 has a classical bulge from the best-fit 
galaxy models, although we cannot completely exclude 
the possibility of this galaxy containing a pseudobulge 
with B/T=0.36. 

The NIR light curves enabled us to isolate the fluxes 
from the host galaxy after ^ 500 days following the dis¬ 
sipation of the X-ray flux. On the other hand, we found 
that there are no significant changes in the light curves of 
the short wavelength bands, supporting the red feature 
of the transient possibly being caused by severe dust ex¬ 
tinction. 

We fit SEDs to the multi-band fluxes of the host galaxy 
which are derived in the light curve analysis. The esti¬ 
mated stellar mass of the host galaxy is log(M*/MQ) = 
9.141 q ®q. The e-folding time scale r is O.IO^oho Gyr 
and the age of stellar population is 0.63tQ43 Gyr. The 
SFR of galaxy is 0.03lQ;g3 Mg/yr. In terms of the sur¬ 
face brightness profile and the stellar mass, this galaxy 
resembles M32, a small companion galaxy of M31. 

We estimated the central Mrh through scaling 
relations. The mass of the SMBH is estimated to be 
10 ® ®'=‘=®-4Mq from Mrh - M*,bui and Mrh - Abui 
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relations for the K band. However, the limit on Mbh 
can be much looser if the host galaxy has a pseu¬ 
dobulge. Future high-resolution, deep imaging should 
be able to unambiguosly distinguish the two possibilities. 
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